Jet quenching parameter q in the stochastic QCD vacuum with 

Landau damping 



Dmitri Antonov, Hans-Jiirgen Pirner 

Institut fur Theoretische Physik, Universitdt Heidelberg, 
Philosophenweg 19, D- 69 120 Heidelberg, Germany 

'66 . 

Q , We argue that the radiative energy loss of a parton traversing the quark-gluon 

o 

CN ■ plasma is determined by Landau damping of soft modes in the plasma. Using this 

U 

CLc idea, we calculate the jet quenching parameter of a gluon. The calculation is done 

in SU(3) quenched QCD within the stochastic vacuum model. At the LHC-relevant 
temperatures, the result depends on the gluon condensate, the vacuum correlation 
i' length, and the gluon Debye mass. Numerically, when the temperature varies from 



o 



X 



T = Tc to T = 900 MeV, the jet quenching parameter rises from q = to approx- 
imately l.SGeVVfm. We compare our results with the predictions of perturbative 



CN ' QCD and other calculations. 

> 

o 
in 



INTRODUCTION 



At RHIC and LHC energies, radiative energy loss is an important mechanism for jet 
energy loss in the quark-gluon plasma It is related to the mean transverse momentum 
(pi) which the parton acquires traversing the plasma. The momentum broadening of the 
■ parton is proportional to the distance it travels, therefore the integrated total energy loss 
from radiation due to the acceleration is proportional to the distance travelled squared Ly 
The energy loss is parametrized by the jet quenching parameter q 

AE = ^CnqLl 

where Cr is the quadratic Casimir operator of the representation R of the parton, and as is 
the strong coupling constant at an appropriate scale. The jet quenching parameter can be 
estimated in a dilute plasma perturbatively: 
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The perturbative temperature dependence of q is determined by the density of scattering 
partners p ocT^ and the differential transport Coulomb cross section, which is cut off in the 
infra-red at the Debye mass. 

A standard nonperturbative calculation of the jet quenching parameter can be done with 
the help of the dipole formalism |2i|. In this formalism, a fake dipole of size L± is constructed 
from the partons in the T-amplitude and in the T*-amplitude, the trajectories of which are 
displaced from each other by the distance L±. The transport parameter {crp^) can then be 
calculated from the cross section of this fake dipole. In the dipole model, for a dipole of size 
L± one has: 

= tr y d^dVih) - V{h + L^)]') = 
= tr y d% (2 - V{h)V\h + Lx) - V\h)Vih + Lx)> . 



Here 



V{h) = V exp 



is the Wilson line of a parton with the impact parametjy: b propagating with the 4- velocity 



= (1, v), = 1, vb = 0, and x{t) = (r, b + vr) [29|. The nonperturbative differential 
cross section to produce a parton with the transverse momentum p± is then given by Fourier 
transform 



da 



dp 



,2 



j rf^L^e'P^^^ tr j d% (2 - V{h)V\h + L^) - V\h)V{h + L^)) 



The expectation values (. . .) of the Wilson lines have to be evaluated with the target ground 
states. The excited states are summed over. The transport parameter can then be obtained 
by differentiation 

= j d'L^_ (-Vie*^"-^) tr j <i=() (2 - V(b)l/t(b + Lx) - V'(b)V(b + Lx)> . 

After a partial integration, one sees that because of the p^-integration only the L^^-dependent 
part, i.e. the first-order term in the dipole cross section is relevant for p^-broadening cal- 
culation. When one takes the expectation value in the medium, the density of scattering 
partners in the medium enters. One can replace the traced products of Wilson lines, tr VV^ 
and tr V^V, by two gauge-invariant Wilson loops, which facilitates the calculation. For 
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example, one has 



tr\/(b)yt(b + Lx) = tr Pexp (^tgv^ J^^ dr [A^(x(r)) - A^(y(r))]| , 

where y(r) = (r, b + — vr). One can approximate two long parallel lines x(r) and ?/(r) 
by a contour, which closes asymptotically at r = ±00, where the interactions vanish. This 
yields the Wilson loop 



tr1/(b)V"^(b + Lx) ~trP exp 
defined at the united contour 



1 

ig I dsZf,{s)Af,{z{s)) 



. xJt(s)), s G fO, 
z^{s) = { ^ V ' 27 ^ ^j^g^g ^^^^ ^ tan(27r(s - 1/4)). 

2/^(r(s)), s G (i,l) 

At high energies, the S'-matrix for dipole scattering is mostly real- valued because of purely 
absorptive interactions (Pomeron exchange). Other contributions have been discussed in the 

fin 

literature [3|, U]. In the calculation, which we will present in this paper, the S'-matrix has 
both real- and imaginary-valued parts, therefore we have to guarantee the reality of the 
cross section by using the real- valued part of the expectation value of the Wilson loop. The 
real- valued part of the expectation value of the Wilson loop in the medium can be related 
to the transport parameter q as follows (see e.g. Ref. [s]): 

{ReW^-l^)=e.p^-^L,Liy (1) 

As a tool to obtain the Wilson-loop expectation value we will use the so-called stochastic 
vacuum model (SVM) {g] in Minkowski space-time. This model has been applied to high 
energy hadron-proton and 7* -proton scattering j?] with two Wilson loops representing pro- 
jectile and target particles. Recently, a single Wilson loop near the light cone has also been 
studied within the SVM in Ref. js]. This calculation is tricky since the expectation value 
of the Wilson loop on the light cone is strictly zero in vacuum, the limiting procedure of 
approaching the light cone, however, can give the quark-antiquark confining interaction in 
the light-cone Hamiltonian. 

In Fig. [1] we show the arrangement which we use for the following calculation. The 
long side of the contour of the Wilson loop is directed along the light cone. The mean 
is related to the transverse size L_i_ of the contour. A counterintuitive property of 
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Figure 1: The contour of the Wilson loop (Ly x L±) of a gluon. 

Eq. ^ is its exponential fall-off, since normally an expectation value of a Wilson loop in 
Minkowski space-time contains the imaginary unity i in the argument of the exponent j^. 
We will demonstrate that the exponential fall-off is a consequence of Landau damping of soft 
modes in the quark-gluon plasma and calculate through this mechanism the jet quenching 
parameter q. 

n 

Recently, a calculation of q at strong coupling has been done for A/" = 4 SYM in Refs. [5||, 
with the help of the AdS/CFT-correspondence. It yields g oc T^, which is a consequence 
of conformal invariance. However, it may be that at LHC-relevant temperatures, the QCD 
conformal anomaly still plays an important role. Lattice results on {e — 3p) ^ indicate that 
the energy density e and the pressure p deviate from pure Stefan-Boltzmann behavior. For 
this reason, one may expect a violation of conformal invariance and an explicit dependence 
of q on the gluon condensate, which arises in the SVM. Of course, the chromo-electric 
gluon condensate vanished at high temperatures after deconfinement, but there can remain 
a chromo-magnetic condensate, and nonvanishing correlations between chromo-magnetic 
field strengths in the quark-gluon plasma. One of our main objectives in this paper is to 
investigate their effect on the jet quenching parameter q. 

The paper is organized as follows. In section II, we review some basics about the SVM at 
T = and T > Tc. We also introduce there an effective local formulation of the SVM, which 
will be technically important for the subsequent analysis. In section III, we incorporate 
Landau damping into the model and calculate the jet quenching parameter. In section IV 
we summarize the main results of the paper. 
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II. GENERALITIES ON WILSON LOOPS AND THE SVM 

We start considering the expectation value of a Wilson loop in Euclidean space-time at 
zero temperature (T = 0). In an arbitrary representation R of SU(Ai'c), it reads (w{C)^ = 
(ti V exp(^ig ^ dXfj, A fj^^ . In the SVM, one uses the non-Abelian Stokes' theorem and the 
cumulant expansion to write the Wilson-loop average as follows 



(w{C)j ~ tr exp 



^ ^ ' da^^{x) I daxp{x') {F^,^{x)Fxp{x')) 



2! 4 

S(C) S(C) 



(2) 



Here, S(C) is the surface encircled by the flat contour C. Averages involving gluonic 
field strengths are calculated in the SVM by using the two-point correlation function 
{F^jj{x)F\p{x')) , which has been fitted to lattice data |9||. To simplify notations, we suppress 
two phase factors between the points x and x' in the fundamental representation. Fixing 
the Fock-Schwinger gauge with the reference point in either of these two points, one can see 
that these phase factors are anyhow equal to 1r when the path between the points x and x' 
is the straight line. The approximate equality in Eq. ([2]) is due to the use of the cumulant 
expansion, which gives the correct area law for the heavy quark-antiquark Wilson loop. It is 
supported by lattice data and states that the amplitude of the two-point irreducible 

gauge-invariant correlation function (cumulant) of F^^?, dominates over the amplitudes of 
higher-order cumulants, which are therefore neglected. The factor 1/2! in Eq. ((21) is due to 
the cumulant expansion, and the factor 1/4 is due to the non-Abelian Stokes' theorem. 
Within the SVM, one parametrizes the nonperturbative part of the two-point cumulant 



as follows [6|, 

+ ("A'^i.p - Up5^x) + {upS^x - uxSpp)] Di{u^) \ . (3) 

Here, D and Di are dimensionless functions of the distance u = x — x' normalized by the 
condition -D(O) = Di(0) = 1. The gluon condensate defines the amplitude of this correlation 
function. The function D is responsible for the confining properties of the non-Abelian gauge 
theory. The function Di describes Abelian-like self-interactions of the Wilson loop, which 
do not lead to an expectation value of the form of Eq. ([T]). High-energy scattering data yield 
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0.74 in the vacuum [l5|. We will take k = 1, i.e. disregard the small contribution of 



the function Di altogether. Lattice simulations [lO|-[l^ yield an exponential fall-off of the 
function D: 



-tj,\u\ 



with the inverse vacuum correlation length 



12| 



fi = 894 MeV. 



One can rewrite the Wilson-loop average using the the surface tensor 



S(C) 



The two-dimensional vector ^ = ^2), where G [0, 1], parametrizes differential elements 
da^i, on the surface and w^{^) points towards each differential surface element. Using 

the correlation function D one obtains the characteristic form of the Wilson-loop average: 

Cr 



W{C) 



exp 



(4) 



48(iV2 _ 

In particular, for a contour C whose temporal extension is much larger than its spatial 
extension, one obtains an area law depending on the string tension in the fundamental 



with Cp 



(5) 



For Nc = 3, we take the standard value of the string tension a = (440 MeV)^ to estimate 
the gluon condensate g'^ {{F^uY) — 3.55 GeV^. 

We see from Eq. ((41) that the SVM essentially suggests a representation of the Wilson-loop 
average in terms of an effective local field theory of the field strength tensor 



with the action 



'S'eucI [F] 



1 



(6) 



(7) 



Here, /C is a well-defined local operator, which reproduces Eq. (Jll) after the Gaussian inte- 
gration in Eq. ([6]): 



iW{C) 



exp 



d%Tjf,^{x)IC ^{x -y)Y:^y{y) 



(8) 
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where /C ^(x) is a Green function of the operator /C: /C/C ^ = S{x). Comparing Eqs. ([4]) 
and dHj), we find that the kernel )C~^ has an exponential fall-off: 

^ ^ ~ 6(iV2 - 1) ^ ■ 

The operator JC itself can then be obtained through Fourier transformation. In the coordinate 
representation, it reads 

-1 ( 92 \ 5/2 

At finite temperature (T 7^ 0), the 0(4) space-time symmetry is broken down to the spa- 
tial 0(3) symmetry, and the correlation function ([3]) splits into the following three functions 

(Bf (x)£)J(x')> , {Et(x)Bl(T!)) , (i?f (x)f?J,(i')> , 



which were simulated on the lattice in Refs. [UJ, ll2|. In the deconfined phase (T > Tc), the 
chromo-electric condensate vanishes, and so does the correlation function (x)£'^(x')). 
Furthermore, the amplitude of the mixed electric-magnetic correlation function is by an order 
of magnitude smaller than the amplitude of the magnetic-magnetic correlation function 



For this reason, we disregard the correlation function (Ef'(x)Bi(x')) with respect to the 



function {B^{x)BI{x')) , which can be parametrized as follows [SOfl: 



{Bt{x)BUx')) = ^^^§^f%kD^{n'), (11) 

D^(^2)^g-MT)H. (12) 

The temperature dependencies of the gluon condensate and of the inverse vacuum correlation 
length will be discussed in the next section. 



III. EVALUATION OF q THROUGH LANDAU DAMPING 

We want to calculate now the Wilson-loop average in the gluon plasma, i.e. in a thermal 
environment. At finite temperature, in the Euclidean space-time, the contour of the loop has 
the orientation shown in Fig.[Tl namely it is oriented at 45° in the (0:4, a;3)-plane. Due to the 
a;4-periodicity, the contour O = Ly x L± effectively splits into pieces, whose extensions along 
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the 3rd and the 4th axes are P = 1/T. Such pieces will be referred to as strips. A point lying 
on the surface of the strip closest to the origin can be parametrized by a vector-function 



t^ = (0,0,l,l), r^= (1,0,0,0), 6,2 6 [0,1]. 



(13) 



Since we have argued to disregard the mixed correlator, only the contribution of the chromo- 
magnetic field to the effective action should be taken into account, which yields [cf. Eq. 1^^] 



SEud[B] = J d^x {B11CB'^_ + iBim 
The Wilson-loop average of one strip reads 
«-stip> = tr 1 VB^.e 



13 j 



(14) 



-■5euci[-B] 



exp 





1 d'^x j 


4 j 





d'^x j rf^?/Si3(a;)/C ^ {x - y)Y.i'j,{y) 



(15) 



Note that, for /C ^ given by Eq. (I9|), this result indeed agrees with the one following from a 
direct application of the cumulant expansion, 



tr exp 



<-^np> = tr ^exp j da^.B^^t'^^ ^ ^ 

dai3{w) I dai^iw'yt'' {B^{w)Bl{w')) 



' 2 



when one uses (-B2 (w)i?2(u'')) in the form of Eq. (fTTll . 

In the full Wilson-loop average (w^^i^^ , also interactions between different strips take 
place besides self-interactions of one strip. Let Xk denote the interaction between two strips 
separated from each other by the distance [3k: 



Xk = ^J dauiw) J dai3{w')JC ^{w-w'), 



where w'^ = w^{C,[,^2) + (0, 0, 0)^. Here, the relative distance between the points, 
|Aw(/c)| = \w — w'\, is 

|Aw;(A;)| = 

1/2 



[m'l - 6) + m + im - 6)] + o [lk^^ - 6)'] \ ^ f3VWT2h^+2^, m 
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In the final form of Eq. (flGll . we disregard O [L']_{^2 — ^2)^], since L± p. If this term had 
been retained, the final result for q would depend on the size of the color dipole, L±. The 
full number of strips (i.e. the maximal value of k) is 

n = (17) 

Therefore, the overall contribution to the Wilson-loop average reads 

n— 1 i n— 1 

i=0 k=0 k=0 

For the transport parameter q we need a calculation of the real- valued part of the adjoint 
Wilson-loop expectation value in Minkowski space-time. We would like to discuss this 
expectation value first, before we make the analytic continuation. In Minkowski space-time, 
in general (W^^'^^{C)'^ = rje^^^ . If one interpretes the Wilson loop as an interaction of 
the fake dipole with the medium, then there can be two medium effects. The fake-meson 
interaction can generate real- and imaginary- valued parts of the S-matrix. Absorptive effects 
[cf. Eq. ([1])], namely the reduction of the intensity of the beam when it traverses the medium, 
are described by rj. In addition, there is real- valued phase shift 5 of the wave-function of the 
projectile, which it acquires propagating through the target. In our calculation, both terms 
will be present, but only 77 will contribute to {(Tp\). That is because the real- valued part of 
the S-matrix oc cos(25) differentiated twice with respect to the dipole size and evaluated 
then at = vanishes. With these preliminaries in mind we can now make the important 
analytic continuation to Minkowski space-time. In Minkowski space-time, the action (fT4l) 
goes over to 

5Mink[S] =1 j d\ {B^ICB^ + B^f'Ei^) , (19) 

and Eq. ([TSll becomes 

(ReW^i^^°^ip) = trRe j pS^e-^^-^I^' = 

= Re exp J (fx j d^yJ:is{x)IC~^{x - y)i:i3{y) ■ (20) 

Until now, this average does not have an exponential fall-off. The exponential fall-off only 
appears when we take into account the scattering partners in the medium as external sources. 
The correlation functions must mediate an interaction with the gluons in the heat bath, 
otherwise there is no absorptive scattering. These gluons can be summed over, see Fig. [2l 
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They will then appear as a gluon polarization insertion into the correlation function of the 
field strengths. The imaginary-valued part of this polarization operator comes from the 
on-shell gluons in the intermediate state, which contribute to the absorptive phase. Since 



the imaginary-va 
gluon plasma 



\ ued part of the gluon self-energy is related to the soft modes of the quark- 
17l |. one sees that the energy loss is essentially linked to Landau damping 
of gluons. The soft background gluons with energies u = \po\ <^ |p|, parametrized by 
the SVM correlation function, can be absorbed by thermal gluons, i.e. gluons from the 
thermal bath. Accordingly, the stochastic QCD vacuum becomes modified in a way, which 
is not accessible in Euclidean simulations on the lattice and, therefore, is not encoded in the 
parametrization of the D^-function, Eq. f |T2l ). The imaginary part of the thermal loop for 



the gluon self-energy reads in the hard-loop approximation [3| 



— ninj5{uo — np) 



' 2|p| 



p2 p2 2 \ p^y v 



(21) 



Here mn = gT^^NjS is the Debye mass of the hard thermal gluon in quenched QCD under 
study. The separation of the momentum scales of the correlation function and the thermal 
gluons is rather subtle. The correlation function contains spatial momenta of order n, but 
very small energies u <ti mo, whereas the thermal heat bath contains gluons with energies 
and spatial momenta of order T. For ultrasoft space-like gluons with u"^ <ti p^, only the 
transverse part of Eq. (i2T]l survives: 

Im Ilij{uj,p) [piPj - p^Sij) ■ V{uj,p), 

The linear component of the field-strength tensor (~ diA"^) satisfies the relation 

Ff-\p)F\r^-^) = 2 {P%-P^P^) A"(P)4(-P)- 

Polyakov [l^ has argued that, in the course of integration over hard modes, the cubic (~ A"^) 
and the quartic (~ A'^) terms appear in such a way that they supplement the square of the 
field-strength tensor to the standard non-Abelian form. For this reason, the spatial part of 
the effective action ([7]) in Minkowski space-time may be modified by Im njj(u;,p): 
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Figure 2: The Wilson-loop surface interacting with the hard gluons (drawn curve) through the 
soft ones (wavy curves). The imaginary part of the soft-gluon polarization operator corresponds to 
Landau damping. For simplicity, the splitting of the surface into strips is not depicted. 



d'x {Ft^ [IC{x) - zV{x)] F^^ + i^^rS,,} . 



(23) 



Accordingly, Eq. (1201) changes as 
(ReiySip> = Reexp 



^^ f d^x f d^yj:^,{x) 



Si3(y) 



xy 



Re exp 



d^yT.i3{x) 



V-tJC 



Sl3(y) 



xy 



(24) 



As was discussed, the term proportional to ilC in Eq. (1241) describes the non-absorptive part 
of the interaction between the fake dipole and the medium. It does not contribute to q 
and will henceforth be disregarded. The remaining real-valued part of Eq. (i24l) describes 
jet quenching. The interaction of the Wilson-loop surface with the hard gluons is depicted 
in Fig. [21 There, the closed line is the hard thermal loop, and the wavy lines are soft 
gluons, which mediate the interactions. The chromo-magnetic condensate and the hard 
gluons represent soft and hard components coexisting in the gluon plasma. 

To handle the complicated Eqs. (I23l) and (1241) . we have approximated P(co', p) by a func- 
tion V{p) of the full 4-momentum in the following way [Sl]- We denote uj = CIpL where 
C <^ 1. Then |p| = \l and we have 



M'^{T) 



<{i-e) 



mliT). 



(25) 
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An important property of the SVM for high-energy scattering in Minkowski space-time 
is the same exponential fall-off with transverse distance of the surface-surface correlations as 
in the Euclidean space-time [cf. Eq. (l20l)]. For this reason, one may calculate the operator 
( yc^+pO xy interest in the Euclidean space-time, where V{p) = ^ and, according to 
Eq. PI), 



[p2 + ;x2(T)]5- 

In this equation, 

'27rV(T)^72((F;j2) 



M(T) 



2 

T 



contains the temperature-dependent gluon condensate ^{F^ivf'^j- ^iid the inverse vacuum 
correlation length /i(T). We therefore obtain 



[ 4 y J J (27r)4 + M4(T)J ^ ^ 

One can now reduce the 4-dimensional momentum integration to a 1-dimensional integration 
by using the formula 



O/dpi) = / - w'Dfip) 



(27r)4 47r2|?i;- w'l ^0 

with the Ji-Bessel function. According to Eq. (fT3ll . the product of two infinitesimal surface 
elements reads 

dar3{w)dai^{w') = {f3L^fd\d\'. 

Furthermore, one can simplify the area-area correlation inside one strip (A; = 0), where the 
distance between the points depends only on the variable x introduced in Eq. (fTBl) : 

\w - w'\ = \Aw{0)\ = (3\x\V2. 

Equation (l26l) yields the following expression: 

V =-l./RplvMink \ CnAf{T)M'{T) dx dpp' J,{p(3\x\^) 

Xo m ^ne w^.,,,,^/ ^^^^^ p ^ J ^ ^ ^ M{T)M\T) ' 

(27) 

The I/^-dependence of this result is consistent with the color transparency of a color-neutral 
dipole of size L±_ interacting with the gluon plasma. Furthermore, the x-integral in Eq. (l27l) 
does not diverge at a; = 0, since this singularity is cancelled by the Bessel function. 
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By using Eqs. ([T|), (fTTI) . and (fTSll . we can now determine q: 

4^2 ■ • 



fc=0 

where 



Xk — 77r^5 P-t^j 



167r2 ^ ^ J.iVk^T2k^T2^ Jo v\v^ + l^\T)f + ^iT)M\Ty 
Here, we have used Eq. ( fT6l) for the distance |Aw(A;)| between two points belonging to the 
strips separated from each other by the interval (3k. Our final result reads 

. _ CrU{T)M^{T) _ dx r dpp^Ji [ppVk^ + 2kx + 2x^) 

^" 47r2n ^^"^ ^ v/A;2 + 2kx + 2x^ Jo p\p^ + li\T)f + M{T)M\Ty 

(28) 

Let us consider a gluon traversing the medium, for which R is the adjoint representation, 
and Cr = Nc = 3. Furthermore, we assume that the critical temperature Tc = 270 MeV 
in SU(3) quenched theory [2^. The temperature behavior of q and the growth of q with T 
depend on the functions ((F^j,)^)^, /i(T), and M(T). To see whether this dependence is 
strong or not, we determine g(T) for two possible choices (I and II) of these functions. 

In case I, the nonperturbative value of the thermal strong coupling gijnn.T) is fixed 



by self-consistency [21]: g{mD,T) = 2.5. This choice defines M{T) through Eq. ( l25ll with 
rriD = g{fnD,T)T. As suggested by the low-temperature lattice data 11, [l2|, the inverse 



vacuum correlation length can be taken constant up to the temperature of dimensional 



reduction [22|, T^ r. — 2Tc = 540 MeV. Above this temperature, all dimensionful quantities 
become proportional to the corresponding power of T. Therefore, /i(T) = const until Td.r., 
whereas fi(T) oc T at T > Td.r.. The proportionality coefficient in the last equation can be 
fixed by the continuity of yu(T) across Td.r.: 

/i(T) = 1.66T at T > Td.r.. 



The temperature dependence of the chromo-magnetic condensate was calculated in Ref. 
9' {{Fp% = g' {{FP') coth f^) at T,<T<Td.r.. 



23| 



V2T/ 

It is also nearly constant up to Td.r.. At T > Td.r., g^ {i.Fp'^)rj. oc T'^, and the proportionality 
coefficient again follows from continuity of g"^ ((T^j^)^)^ across Td.r.: 

/((0'>t = 61-46T4 at T > Td.r.. 
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Figure 3: The chromo-magnetic condensate as a function of temperature at Tc < T < 900 MeV in 
cases I and II. 

In case II, we define the parameter M(T) through the perturbative 1-loop strong cou- 



pling 



22|: 



g~\T)=2bolnj, where &o = ^ 



11 

Nc=3 IGvr^ 



A ^ 0.104T, 



Furthermore, high-temperature lattice data 



22| suggest the parametrization 



/i(r) = i.04/(r)T, 

where the coefficient has been fixed by the value /i(Tc) = fi. The same lattice data suggest 
also a temperature-dependent spatial string tension 

a{T) = [cg^{T)T]'^, where c = 0.566. 

Note that such a value of the coefficient c corresponds to a spatial string tension 
cr(Tc) = (485 MeV) ^, that is slightly larger than the standard zero-temperature value 
cr(0) = (440 MeV)^. The temperature-dependent chromo-magnetic condensate can then 
be obtained from Eq. ^ = 3: g^{iFp^)^ = ^iJ.'^{T)a{T). For illustration, we 

present the chromo-magnetic condensates for both cases I and II in Fig. [3l 

Because of the much smaller density of scattering partners in the hadronic phase, typical 
values of g at T < Tc are ~ 0.01 GeV^/fm, which is by two orders of magnitude smaller than 
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in the quark-gluon plasma [2j|. For this reason, we set q(T) = q(T) — q(Tc), i.e. q(Tc) = 0. 
Furthermore, we truncate n = 10 ^ 1 and use a small value C = 0.1 ^ 1. We then obtain q 
as a function of temperature in the interval 270 MeV < T < 900 MeV for both cases I and 
II, cf. Fig. m Both calculations of q follow a temperature dependence oc T^. We find 

q{T)cLi = 0.16(T/T,)3 GeVVfm, q{T)cLu = 0.2Q{T/T,)' GeV'/fm. 



These fits are inspired by Eq. f l28ll . Indeed, counting the powers of T in this equation at 
high temperatures, one sees that 

q(xT' at T>Td,.. 

Furthermore, at T = 900 MeV we have tested the sensitivity of our results against an increase 
of the number of strips n. An increase of the number of strips leads to a small increase of 
g(T), namely 

g(900 MeV)i"=i°^ = 1.26 GeVVfm, g(900 MeV)^";;i°^ = 1.39 GeVVfm; 

g(900 MeV)^"=if = 1.78 GeVVfm, g(900 MeV)i"=if = 1.98 GeVVfm. 

We have also addressed the dependence of q{T) on the parameter (, which measures the 
smallness of the energy of thermal gluons with respect to their characteristic spatial mo- 
menta. For small values of ( under consideration, specifically ^ < ~ 0.58, the function 
M^(T), Eq. (l25ll . describing the imaginary part of the hard thermal loop polarization op- 
erator, is a monotonically increasing function of (. Therefore, since M'^{T) enters both the 
numerator and the denominator of the resulting Eq. (i28ll , only a numerical analysis can say 
whether q{T) increases or decreases with the increase of (. We have depicted the results of 
such an analysis in Fig. [5] for case I and in Fig. [6] for case II. In case I, we observe a slow 
monotonic decrease of g(T) with the increase of (. In case II, we observe a small increase 
of g(T) when ( varies from 0.1 to 0.3 and an essentially constant behavior of q(T) with the 
further increase of (. In both cases, the dependences of q(T) on ( are very weak, which 
indicates that the predictions of our model are rather stable. 



IV. CONCLUDING REMARKS 



In this paper, we have evaluated the jet quenching parameter q in SU(3) YM theory. Our 
calculation is based on a nonperturbative approach, represented by the stochastic vacuum 
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Figure 4: The jet quenching parameter q{T) in case I (lower full curve) and II (upper dashed curve), 
and the interpolating fitted curves ~ T3 in both cases. 
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Figure 5: The jet quenching parameter q{T) in case I for various values of the parameter (. 
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Figure 6: The jet quenching parameter q{T) in case II for various values of the parameter C- 



model, which describes the scattering of a fast parton off a hard thermal gluon in the leading 
approximation oc {g'^ {{F^luY) rpY ■ In our scenario of jet quenching, the gluon plasma has two 
components. The chromo-magnetic condensate describes the soft component of the gluon 
plasma, which has been called "epoxy" in the works of the Stony Brook group 2^- The hard 
component of the plasma (with momenta larger than the inverse vacuum correlation length) 
is represented in our approach by the hard thermal loop effective theory. In some sense, such 
a picture of a two-compone nt g luon plasma resembles the two-component Landau model of 



superfluid quantum liquids 



In our model, jet quenching originates from Landau damping of soft gluons by the on- 
shell hard thermal gluons. The analytic result for q(T) is given by Eq. (i28l ). Numerically, we 
have calculated q{T) for two alternative parametrizations of the chromo-magnetic condensate 
and the vacuum correlation length, thereby testing the stability of the result against possible 
uncertainties of these quantities. The numerical results are plotted in Fig. [4] along with the 
corresponding fitting curves. Our values of q are somewhat larger than the values obtained 
in pQCD, qp o n n = 1-1 -r- 1.4 GeV^/fm jl| and closer to those of other recent nonperturbative 



calculations 
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27|, where g = 1.0 ^ 1.9GeV^/fm. However, due to the difference of our 



model from that of Ref. 



27||, our result Eq. (i28l ) differs parametrically from the result of 
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that paper. 

Note finally that our values of the jet quenching parameter are close also to the effective 
values in A/" = 4 SYM corresponding to QCD. Indeed, analyzing possible matchings of 
the two theories, Gubser [2^ has argued that the result of Ref. would correspond to the 
values of the jet quenching parameter of a heavy quark in QCD q = 0.61 ^ 2.6GeV^/fm. 
Since the result of Ref. is not proportional to the Casimir operator Cr, the values of the 
jet quenching parameter for a gluon should be the same. These values are of the same order 
of magnitude as those obtained in the present paper. 

Acknowledgments 

We are grateful to D. D. Dietrich for a collaboration at the early stage of the work. The 
work of D.A. has been supported through the contract MEIF-CT-2005-024196. 



[1] R. Baler, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 483, 291 
(1997); for reviews see: R. Baier, D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50, 
37 (2000); M. Gyulassy, I. Vitev, X. N. Wang and B. W. Zhang, "Jet quenching and radiative 
energy loss in dense nuclear matter," |arXiv:nucl-th/0302077 

[2] J. Dolejsi, J. Hiifner and B. Z. Kopeliovich, Phys. Lett. B 312, 235 (1993). 

[3] O. Nachtmann, Annals Phys. 209 (1991) 436. 

[4] For a review see: S. Donnachie, H. G. Dosch, O. Nachtmann and P. Landshoff, Pomeron 
physics and QCD, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 19, 1 (2002). 

[5] H. Liu, K. Rajagopal and U. A. Wiedemann, Phys. Rev. Lett. 97, 182301 (2006); JHEP 0703, 
066 (2007). 

[6] H. G. Dosch, Phys. Lett. B 190, 177 (1987); H. G. Dosch and Yu. A. Simonov, Phys. Lett. 

B 205, 339 (1988); for a review see: A. Di Giacomo, H. G. Dosch, V. L Shevchenko and 

Yu. A. Simonov, Phys. Rept. 372, 319 (2002). 
[7] A. L Shoshi, F. D. Steffen and H. J. Pirner, Nucl. Phys. A 709, 131 (2002). 
[8] H. J. Pirner and N. Nurpeissov, Phys. Lett. B 595, 379 (2004). 
[9] E. Meggiolaro, Phys. Lett. B 451, 414 (1999). 



19 



[10 

[11 
[12 

[13 



A. Di Giacomo and H. Panagopoulos, Phys. Lett. B 285, 133 (1992). 

A. Di Giacomo, E. Meggiolaro and H. Panagopoulos, Nucl. Phys. B 483, 371 (1997). 

M. D'Elia, A. Di Giacomo and E. Meggiolaro, Phys. Rev. D 67, 114504 (2003). 



For reviews see: A. Di Giacomo, "Non perturbative QCD," arXiv:hep-lat/9912016 ; Czech. J. 



Phys. 51, B9 (2001); "QCD vacuum and confinement," |arXiv:hep-lat/0204001 



[14] G. S. Bah, N. Brambiha and A. Vairo, Phys. Lett. B 421, 265 (1998). 

[15] A. L Shoshi, F. D. Steffen, H. G. Dosch and H. J. Pirner, Phys. Rev. D 68, 074004 (2003). 

[16] D. Antonov, D. Ebert and Yu. A. Simonov, Mod. Phys. Lett. A 11, 1905 (1996); D. Antonov, 
Surveys High Energ. Phys. 14, 265 (2000), section 2.1. 

[17] For a review see: J. P. Blaizot and E. lancu, Phys. Rept. 359, 355 (2002). 

[18] V. V. Klimov, Sov. J. Nucl. Phys. 33, 934 (1981); Sov. Phys. JETP 55, 199 (1982); H. A. Wel- 
don, Phys. Rev. D 26, 1394 (1982). 

[19] A. M. Polyakov, Gauge Fields and Strings (Harwood Academic Publishers, Chur, 1987), chap- 
ter 2.3; see also: D. Antonov and H. J. Pirner, Eur. Phys. J. C 51, 633 (2007). 

[20] For a recent review see: P. Petreczky, Eur. Phys. J. C 43, 51 (2005). 

[21] J. Braun and H. J. Pirner, Phys. Rev. D 75, 054031 (2007); D. Antonov, S. Domdey and 

H. J. Pirner, Nucl. Phys. A 789, 357 (2007). 
[22] G. S. Bah, J. Fingberg, U. M. Heller, F. Karsch and K. Schilling, Phys. Rev. Lett. 71, 3059 

(1993); F. Karsch, E. Larmann and M. Liitgemeier, Phys. Lett. B 346, 94 (1995); G. Boyd, 

J. Engels, F. Karsch, E. Larmann, C. Legeland, M. Liitgemeier and B. Petersson, Nucl. Phys. 

B 469, 419 (1996). 
[23] N. O. Agasian, Phys. Lett. B 562, 257 (2003). 

[24] For a recent reference see: J. Casalderrey-Solana and X. N. Wang, Phys. Rev. C 77, 024902 
(2008). 

[25] G. E. Brown, C. H. Lee, M. Rho and E. Shuryak, Nucl. Phys. A 740, 171 (2004); G. E. Brown, 
L. Grandchamp, C. H. Lee and M. Rho, Phys. Rept. 391, 353 (2004). 

[26] See e.g.: E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part II (Pergamon, Oxford, 
1980). 

[27] A. Majumder, B. Miiller and X. N. Wang, Phys. Rev. Lett. 99, 192301 (2007). 
[28] S. S. Gubser, Nucl. Phys. B 790, 175 (2008). 

[29] Throughout the paper, we use the notation = A'^t'^, where t'^'s are the generators of the 



20 



group SU(A^c) in a given representation R. The trace "tr" is normalized by dividing over the 
trace of the unit matrix in the representation R. Therefore, with this definition, tr l/j = 1. 
The normahzation of the function can be proven by taking the (T — > 0)-hmit and imposing 
for the correlation function (^E^{x)E\.{x')^ the same parametrization (fTTI) with some other 
function instead of . Due to the unbroken 0(4)-invariance at T = 0, one has Z?^(0) = 
D^(0). On the other hand, using the parametrization (fTTI) . one has 

i=l 

This proves that D^{Q) = 1. 

For brevity, we use the same notations for a function or an operator and its Fourier image. 
What is implied in each particular case is clearly defined by the corresponding argument. 



